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Abstract 

We propose a macroscopic modeling framework for a network of roads and multi-commodity traf¬ 
fic. The proposed framework is based on the Lighthill-Whitham-Richards kinematic wave theory; more 
precisely, on its discretization, the Cell Transmission Model (CTM), adapted for networks and multi- 
commodity traffic. The resulting model is called the Link-Node CTM (LNCTM). 

In the LNCTM, we use the fundamental diagram of an “inverse lambda” shape that allows modeling 
of the capacity drop and the hysteresis behavior of the traffic state in a link that goes from free flow to 
congestion and back. 

A model of the node with multiple input and multiple output links accepting multi-commodity traffic 
is a cornerstone of the LNCTM. We present the multi-input-multi-output (MIMO) node model for multi- 
commodity traffic that supersedes previously developed node models. The analysis and comparison with 
previous node models are provided. 

Sometimes, certain traffic commodities may choose between multiple output links in a node based 
on the current traffic state of the node’s input and output links. For such situations, we propose a 
local traffic assignment algorithm that computes how incoming traffic of a certain commodity should be 
distributed between output links, if this information is not known a priori. 


Keywords: macroscopic first order traffic model, first order node model, multi-commodity traffic, dynamic 
traffic assignment, congestion games, managed lanes 


1 Introduction 


Traffic simulation models are important tools for traffic engineers and practitioners. As in other disciplines 
such as climatology, population dynamics, etc., traffic models have helped to deepen our understanding of 
traffic behavior. They are widely used in transportation planning projects in which capital investments 
must be justified with simulation-based studies (Caltrans, 2015). Recently, with the increased interest in 
Integrated Corridor Management (ICM) and Decision Support Systems (DSS), traffic models have found a 
new role in real-time operations. The requirements for models used in the real-time context, both in terms 
of execution speed and modeling accuracy, are significantly higher than in the planning context. A model 
used within an on-line traffic prediction framework, for example, must adapt to incidents, weather changes, 
special events, as well as to the normal day-to-day variations in demand. It must also be capable of capturing 
all aspects of the transportation system that are relevant to the operator’s decision making process. These 
may include arterials as well as freeways, HOV and HOT lanes, ramp meters, variable speed limits, etc. 
They may also include multiple modes of travel (buses, trains, cars) and their respective emission signatures. 
Traffic models are also a central component of state estimation (Wright and Horowitz, 2015), in which an 
entire ensemble of simulations covering a range of possible states must be executed in real time. 
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The ICM and traffic DSS systems that have been deployed in the U.S. are all driven by microsimulation 
models. These models simulate the interactions between individual vehicles, and are therefore much more 
computationally expensive than macroscopic models such as the Cell Transmission Model (Daganzo, 1994, 
1995). One of the reasons for this discrepancy is that macroscopic models have lagged behind their com¬ 
mercial microscopic counterparts in their ability to simulate complex configurations such as HOV lanes, 
arterial/freeway interactions, and multi-commodity flows. 

This paper attempts to narrow the gap by introducing some new features for macroscopic models. These 
features are gathered in a model we call the Link-Node Cell Transmission Model, which, as its name implies, 
is a generalization of Daganzo’s original CTM. Next we describe each of the new features. 

First we introduce a link model that incorporates the notions of capacity drop and hysteresis that have 
been observed in empirical studies (Srivastava and Gerolimini, 2013). The term “capacity drop” refers to the 
observation that the maximum flow that can be measured directly downstream of a bottleneck is lower if the 
bottleneck is active than if it is not (Coifman and Kim, 2010). Hysteresis in traffic dynamics refers to the 
fact that recovery from congestion usually follows a different path in the flow/density plane than the descent 
into congestion (Saberi and Mahmassani, 2013). Both of these together can be modeled with a special form 
of the fundamental diagram commonly referred to as the “inverted lambda” form (Koshi et al., 1983) 1 . Of 
note is that the “inverted lambda” destroys the static functional relationship between flow and density, and 
therefore requires an extra state, which we will call the “congestion metastate”. 

The second contribution of the paper relates to the treatment of network nodes. Node models determine 
how congestion will travel and distribute from one link to upstream links. For example, the node model 
determines whether congestion on a freeway will remain on the mainline or spill through the onramps onto 
the streets. We present a new node model which unifies the existing models of Tampere et al. (2011) and 
Bliemer (2007). This model uses input link priority parameters and can be applied to multi-commodity 
traffic. Input link priorities define how output supply should be allocated for incoming flows. In Tampere 
et al. (2011), input link priorities were assigned implicitly. Our node model admits arbitrary input link 
priorities. 

Third, we propose a relaxation of the first-in-first-out principle (FIFO) that is common to most if not all 
macroscopic models. The FIFO principle states that all vehicles in a link must travel at the same speed, 
irrespective of their type or destination. This is not a good approximation in the case of multi-lane traffic, 
in which each lane is allowed to move at its own speed. Applied to a freeway, it means that congestion 
backing into the mainline from an offramp will immediately block all lanes. In reality vehicles headed for a 
congested offramp will tend to queue in the slow lane, while leaving the remaining lanes to continuing traffic. 
The solution to this problem adopted by authors such as Shiomi et al. (2015) is to assign a state variable 
to every lane in every link. We propose here a simpler, more computationally efficient scheme in which split 
ratios are adjusted depending on the downstream availability of space. We allow FIFO to be relaxed across 
destinations, while preserving FIFO across vehicle types. 

The fourth contribution of this article relates to traffic behavior at complex junctions such as access points to 
managed lanes (e.g., HOV and HOT lanes). Drivers are allowed to choose whether or not to use a managed 
lane, and they make that choice based on the current observed state of congestion. For example, the incentive 
to use an HOT lane is greater when the travel time difference between the HOT and the general purpose 
lanes is large. Thus, a driver’s decision to pay for HOT access will depend on the expected travel times 
and the price. Split ratios in this situation are not known beforehand - they are endogenous to the model. 
We develop a split ratio assignment algorithm 2 for this problem in which the split ratios across an HOV 
or HOT gate are computed by the simulator. This same algorithm may be adapted to other situations in 
which turning decisions depend on local conditions, as is the case when arterial traffic routes itself around 
an obstruction. 

1 Other names for this shape used by various authors include “reverse lambda” and “backwards lambda.” 

2 In this paper, the terms “local [dynamic] traffic assignment” and “split ratio assignment” are used interchangeably. 
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The paper is organized as follows. Section 2, describes the Link-Node Cell Transmission Model (LNCTM). 
The main feature of this Section is the “inverse lambda’-shaped fundamental diagram. Section 3 presents our 
new traffic node model, starting with the analysis of the simpler multiple-input-single-output (MISO) case, 
then moving to the general multiple-input-multiple-output (MIMO) case. Then, we introduce the relaxed 
FIFO rule. The proposed node model is used for computing cross-flows into and out of HOV lanes along 
the freeway. Section 4 describes the split ratio assignment algorithm. Section 5 concludes the paper by 
summarizing the results. For the convenience of the reader, the notation used in this paper is summarized 
in the Appendix. 


2 Link-Node Cell Transmission Model 


The LNCTM models traffic flow in a road network consisting of links C and nodes A f, where links represent 
stretches of roads, and nodes represent junctions that connect links. A node always has at least one input 
and at least one output link. A link is called ordinary if it has both begin and end nodes. A link with no 
begin node is called origin, and a link with no end node is called destination. Origins are links through 
which vehicles enter the system, and destinations are links that let vehicles out. 



Figure 1: Fundamental diagram. 


Each link l £ C is characterized by its length and the time dependent fundamental diagram, a flow-density 
relationship presented in Figure 1. A fundamental diagram may be time-varying and is defined by 4 values: 
capacity F;(t), free flow speed vf (t), congestion wave speed un(t) and the jam density nf (t). In this paper 
we assume that densities, flows and speeds are normalized by link lengths and discretization time step; 3 
and that free flow speed vf ( t ) and congestion wave speed Wi(t) satisfy the Courant-Friedrichs-Lewy (CFL) 
condition (Courant et al., 1928): 0 < vf (t),wi(t) < 1. 4 5 6 The values nff {t) = an d nf(t) = are 

called low and high critical density respectively. Unless nT (t) = nf(t), when it assumes triangular shape, 
the fundamental diagram is not a function of density: n/(t) £ (nff {t),nf (t)] admits two possible flow values. 


Each node v £ Af with Al input and N output links is characterized by time dependent mutual restriction 
coefficients {i)jj'{t)Y, 5 input link priorities {pi(t)} and partially defined split ratios {/3fj(t)}, 6 where C is 
the number of vehicle types; i = 1,..., M, j,j' = l,...,N and c = 1 ,... ,C. 

3 Civen original (not normalized) capacity F) specified in vehicles per hour (vph), free flow speed vf and congestion wave 
speed u>i specified in miles per hour (mph), and jam density nf specified in vehicles per mile (vpm), as well as link length A;/:/ 
and discretization time step At, normalized values are F\ = if At specified in vehicles per time period At, vf = vf -fff~ and 
Wj = wi Yf , both unitless, and rtf = nf A xi specified in vehicles. 

4 The CFL condition is the necessary condition for convergence while solving hyperbolic PDEs numerically. 

5 Mutual restriction coefficients determine how flow restrictions in output links influence each other, when the first-in-first-out 
(FIFO) condition us relaxed. This concept is explained in Section 3.4. 

6 Split ratios may also be fully defined or fully undefined. 
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The state of the system at time t is described by the number of vehicles per commodity in each link: 
ni(t) = [■ nj(t ),..., np(t)] , where nf(t) represents the number of vehicles of type c in link l at time t. The 
state update equation for link l £ £ is: 


n l (t + l)=n l (t)+ (frit)-frit)) 


( 2 . 1 ) 


C,in 


where /*»(i) = 


iT 


is the vector of commodity flows coming into link l during this time step, 
is the vector of commodity flows leaving link l during this time step. 


For ordinary and destination links, f[ n (t) is obtained from the begin node: given a begin node with M input 
links, 


M 

/r <n (t) = £/«(*), c = 1,..., C. (2.2) 

i— 1 

For origin links, 

fr n (t)=dt(t ), ( 2 . 3 ) 

where df(t) denotes commodity demand at time t, which is an exogenous input to the model, specified in 
vehicles per discretization step At. 

For ordinary and origin links, f° ut (t) is obtained from the end node: given an end node with N output links, 


N 

/r Wt (*) =£/«(*)’ c=l,...,C. (2.4) 

i=i 

For destination links, 

f[’ out (t) = v{ (fK(t) min 11, c F [ e> 1 , c = 1,..., C. (2.5) 

t z ^ c '=i v i v) n i v ) J 

The values f-j(t) and /(( (t) are computed by the node model that is addressed in detail in Section 3. 7 
For each link l £ £ we will also define a congestion metastate: 

( 0 ni(t)<nf(t), 

Oi{t) = < 1 ni{t) > n/"(t), (2.6) 

[ Oi(t- 1) nf (■ t ) < ni(t) < np(t), 

where ni(t) = 5P C=1 n°(t). This metastate helps determining which constraint of the fundamental diagram 
is activated when we compute the receive function for a link. 

Now we can formally describe the LNCTM that runs for T time steps. 

1. Initialize: 

<(0) := n? i0 ; 

9i(0) := 9i t o; 

t := 0 

for all l £ £ and c = 1,..., C, where n^ 0 and 6i 0 are initial conditions. 

7 More precisely, by the MIMO algorithm with relaxed FIFO condition described in Section 3.5. 
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2. Apply all the control functions that modify system parameters (fundamental diagrams, input priorities) 
and/or system state. A control function may represent ramp metering, variable speed limit, managed 
lane policy, etc. Control functions may be open-loop (if they depend only on time) and closed-loop (if 
they depend on time and system state). This step is optional. 

3. For each link l £ C and commodity c = 1,..., C define the send function: 


sm = 


v{ m i R {] 

d c { ( t ) min j ] 


Fi(t) 


r WE c c 


mm 

Ef: ■ . 


rw 


}■ 


l is an ordinary link or a destination, 
l is an origin. 


(2.7) 


4. For each link l £ C define the receive function: 


Ri{t) 


(1 - 6i{t)) Fi(t) + Oi{t)wi(t) (nf (t) - Yfc= i > 
00 , 


l is an ordinary link or a destination, 
l is an origin. 

( 2 . 8 ) 


5. For each node v £ Af that has undefined split ratios, given its input link priorities { Pi(t)}, send 

functions Sf(t) and receive functions Rj(t ), compute the undefined split ratios according to 

the algorithm from Section 4. 

6. For each node v £ A f, given its mutual restriction coefficients { r //y(i)}, input link priorities {pi(t)} 
and split ratios {/3()(t)}, send functions Sf(t') and receive functions Rj(t ), compute input-output flows 
f£ according to the algorithm from Section 3.5. 

7. For each link l £ £, compute /™(t) using expressions (2.2)-(2.3) and f° ut (t) using expressions (2.4)- 
(2.5). 

8. For each link l £ £, update the state ni(t + 1) according to the conservation equation (2.1), and the 
metastate 9i(t + 1) according to its definition (2.6). 

9. If t = T, then stop, otherwise set t := t + 1 and return to step 2. 


Traffic speed for link l is computed as a ratio of total flow leaving this link to the total number of vehicles 
in this link: 


Vl[t) 


E/Li ir' 

Ec=i "f 

v{ (t), 


if Ef= i > 0: 

otherwise. 


(2.9) 


Defined this way, Vi(t) £ [0,vf (t)]. 


3 Node Model 


As mentioned in Section 2, at each timestep of the LNCTM, traffic flows between links are calculated at the 
node level, such that flows through the node are functions of the state of all links joined at the node (step 5 
of the LNCTM algorithm, see Section 2). For nodes with simple link configurations, such as the single-chain 
case with M = 1,N = 1, these calculations are straightforward, but the present problem requires a node 
model formulation for computing input-output for general values of M, N , and C. 

In Tampere et al. (2011), the authors proposed a list of eight requirements for any general node model. With 
our amendments, this list looks as follows. 
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1. Applicability to general numbers of input links M and output links N. We extend this requirement to 
include the general number of traffic commodities C. 

2. Maximization of the total flow coming out of the node. According to Tampere et al. (2011), it means 
that “each flow should be actively restricted by one of the constraints, otherwise it would increase until 
it hits some constraint”. When a node model is formulated as a constrained optimization problem, its 
solution will automatically satisfy this requirement. However, what this requirement really means is 
that constraints should be stated correctly and not be overly simplified and, thus, overly restrictive for 
the sake of convenient problem formulation. 

3. Non-negativity of all input-output flows. 

4. Flow conservation: total flow entering the node must be equal to the total flow exiting. This require¬ 
ment is automatically satisfied, since we deal directly with input-output flows inside the node. 

5. Satisfaction of demand and supply constraints. 

6. Satisfaction of the first-in-first-out (FIFO) constraint: if a single destination j for a given i is not able 
to accept all demand from i to j, all flows from i are also constrained. 

We believe that in some situations, the FIFO constraint may be too restrictive. It should not be 
completely eliminated, however, but must be relaxed through a parametrization. This FIFO relaxation 
is discribed in Sections 3.4 and 3.5. Parameters that we refer to as mutual restriction coefficients allow 
configuring each pair of node’s output links to anything between no FIFO at all and fully enforced 
FIFO. 

7. Satisfaction of the invariance principle. If the flow from some input link i is restricted by the available 
output supply, this input link enters a congested regime. This creates a queue in this input link and 
causes its demand S) to jump to capacity F t in an infinitesimal time, and therefore, a node model 
should yield solutions that are invariant to replacing 5) with C) when flow from input link i is supply 
constrained (Lebacque and Khoshyaran, 2005). 

Although the invariance principle is especially important in developing numerical methods for contin¬ 
uous time solutions, whereas our traffic model is discrete in space and time, and its time step At is 
not infinitesimal, we retain it in the list of node model requirements and will revisit it in Section 3.3. 

8. Satisfaction of the supply constraint interaction rule (SCIR), which is supposed to represent the ag¬ 
gregate driver behavior at congested nodes. Following Gentile et al. (2007), in Tampere et al. (2011) 
it was proposed to allocate supply for incoming flows proportionally to input link capacities. 

In this paper we propose the concept of input link priorities that can be represented by arbitrary 
nonnegative values (or functions) as SCIR. Input link priorities will define the allocation of the output 
supply for incoming flows. It must be noted that the wrong choice of input link priorities may lead to 
violation of the invariance principle. 


We will add the 9th requirement concerning multi-commodity nature of modeled traffic to this list. 


9. Supply restriction on a flow from any given input link is imposed on commodity components of this 
flow proportionally to their presence in this link. 


As opposed to the first eight, we do not view requirement 9 as a general principle. Rather, we list it here as 
a rule that we follow throughout this paper. 

Further we will provide the mathematical formulation of the listed requirements and build the node model 
that satisfies them. This Section is organized as follows. We start by describing the merge problem — the 
multi-input-multi-output (MISO) node in Section 3.1 to explain the concept of input link priorities. From 
there we move to the general multi-input-multi-output (MIMO) node in Section 3.2. Then, we discuss the 
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relationship of the proposed node model with the node model of Tampere et al. (2011) and also compare 
it with the node model presented in Bliemer (2007) in Section 3.3. In Section 3.4 we propose a way of 
relaxing the FIFO condition and explain it in the case of single-input-multi-output (SIMO) node. Finally, 
in Section 3.5 we proceed to the general MIMO node with relaxed FIFO condition. 


3.1 Multiple-Input-Single-Output (MISO) Node 


First, we consider a node with M input links and 1 output link. The number of vehicles of type c that input 
link i wants to send is Sf. The flow entering from input link i has priority pi > 0. Here i = 1,M and 
c = 1,... ,C. The output link can receive R vehicles. 


Casting the computation of input-output flows ff as a mathematical programming problem, we arrive at: 


max 


( m c 

EE/- 


i 1 


y i—1 C— 1 


subject to: 


(3.1) 


fii> 0, i = 1,..., M, c = 1,..., C — non-negativity constraint; 
fii < Sf, i = 1,..., M, c = 1,.. ., C — demand constraint; 


M 


EE fa< R — supply constraint; 


i—1 c—1 

fti 


Sf 


i = 1, •. •, M, c = 1, 


, C, — proportionality constraint 


E C rc' CC' 

c'—l J il Z^c' = l 

for commodity flows; 

(a) Pi" Ef=i fti = Pi 1 Ef=i ft "l c - such th at < Sf,, /P, a < 

(b) If Ec= i fa < Ec=i s °a then E^=i fti > 


priority constraint. 


(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 


Generally, constraint (3.5) is optional. In this paper, however, we assume that all constraints on input- 
output flows are imposed on commodity flows proportionally to commodity contributions to total demand. 
This constraint can be interpreted as a FIFO rule for multi-commodity traffic within a given link; not to be 
confused with the FIFO rule for multiple output links. 

Let us discuss priority constraint (3.6) in more detail. It should be interpreted as follows. Input links 
i = 1,..., M, fall into two categories: (1) those whose flow is restricted by the allocated supply of the output 
link; and (2) those whose demand is satisfied by the supply allocated for them in the output link. Priorities 
define how supply in the output link is allocated for input flows. Condition (3.6) (a) says that flows from 
input links of category 1 are allocated proportionally to their priorities. Condition (3.6) (b) ensures that 
input flows of category 2 do not take up more output supply than was allocated for them in cases where 
category 1 is non-empty. The inequality in Condition (3.6) (b) becomes an equality when category 2 is empty. 

There is a special case, when some input link priorities are equal to 0. If there exists input link i, such that 
Pi = 0, while fn > 0, then, due to condition (3.6)(a), all input links with non-zero priorities are in category 
2. Thus, if category 1 contains only input links with zero priorities, one should evaluate condition (3.6)(a) 
with arbitrary positive, but equal , priorities: p^ = pi" > 0. 
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If the priorities pi are proportional to send functions Y^=i i = 1,... , M, 8 then condition (3.6) can be 
written as an equality constraint: 


E C fc 

C— 1 J 11 

E?=i-s? 


E C fc 

C— 1 ■'il 

E C cc 

c=l 


E C rc 
c— 1 J Ml 

E C QC ' 
C— 1 °M 


and the optimization problem (3.1)-(3.6) turns into a linear program (LP). 


(3.7) 


For arbitrary priorities with M = 2, condition (3.6) becomes: 



To give a hint how more complicated constraint (3.6) becomes as M increases, let us write it out for M = 3: 


S* s H (*-£*) ^ (*-£*) > 1 (310) 


<3 ' n) 


c 

E re ^ I P3 n P3 

f 31 < max < - R , 

C— 1 


C 


Eti Pi Pl+P3 




P3 


C 


C 


, ( 3 - 12 ) 

c—1 / w \ c—1 / *=1,2 c=l 


As we can see, right hand sides of inequalities (3.10)-(3.12) contain known quantities, and so for arbitrary 
priorities, problem (3.1)-(3.6) is also an LP. For general M, however, building constraint (3.6) requires a 
somewhat involved algorithm. Instead, we present the algorithm for computing input-output flows /() that 
solves the maximization problem (3.1)-(3.6). 


1. Initialize: 

R{ 0) := R\ 

U( 0) := {1,..., M} ; 
k := 0. 


U(k) is the set of unprocessed input links: input links whose input-output flows have not been assigned 
yet. 


2. Check that at least one of the unprocessed input links has nonzero priority, otherwise, assign equal 
positive priorities to all the unprocessed input links: 


Pi ( k ) = 


Pi, 

WW\-' 


if there exists i' € U{k) : pv > 0, 
otherwise, 


where \U(k)\ denotes the number of elements in set U(k). 


Here we assume that Sf > 0 , i = 1 ,, M. 




3. 


Define the set of input links that want to send fewer vehicles than their allocated supply and whose 
flows are still undetermined: 


U(k) 


c 

GU(k) : <Pi(fc) 

C— 1 


m \ 

J2i'eu(k) Pi'(k) j 


• If U ( k ) ^ 0, assign: 


fa 

R(k + 1) 

U(k + 1) 
• Else, assign: 

fl i 

U(k + 1) 


5 8 c , zel7(fc); 

iSE/ C=1 

U(k)\U{k). 




Pi(k) 


R(k) 


'Ei'eu(k)Pi'( k 'l Ec=i 'S'f ’ 


i G C/ (At) ; 


4. If U(k + 1) = 0, then stop. 

5. Set k := k + 1, and return to step 2. 


This algorithm finishes after no more than M iterations, and in the special case of M = 2 it reduces to 


/u 

f°2 1 


max 


min < Si, Si- 


/ pi 

1 P1+P2 


R, R- 


ES =1 s/} 


E^=i E' 


max 


min < ^ 2 , ^2 


/ P2 
\ P1+P2 


I?, i?- 


ES=X}' 


Ec'=l S 2 


(3.13) 

(3.14) 


The following theorem can be trivially proved, and we leave the proof to the reader. Further, in Section 3.5, 
it will become clear that this theorem is a special case of a more general statement. 

Theorem 3.1. The input-output flow computation algorithm constructs the unique solution of the maxi¬ 
mization problem (3.1)-(3.6). 


Example. Suppose we have three inputs and one output, with parameters: 


C = 1 R = 1000 
S\ = 400 pi = i/3 
s 2 = 500 p 2 = 2 /3 
S 3 = 200 p 3 = 0 
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Our solution algorithm for finding the resulting flows proceeds as follows: 

fc = 0 : 


1. 

17(0) = {1,2,3}; i?(0) = 1000; 


2. 

Pi(0) = i/3, pa(0) = 2 /s, P3(0) = 0 


3. 

1 : 

2. 

17(0) = {2}; f 2 = 500 ; i?(l) = 500; 

17(1) = {1,3} 

pi(l) = V3, p 3 (l) = 0 


3. 

o . 

U( 1) = {1}; fi = 400 ; R( 2) = 100; 

17(2) = {3} 

z . 

2. 

p 3 (2) = 1 


3. 

#(2) = 0; / 3 = 100 ; 17(3) = 0 



3.2 Multiple-Input-Multiple-Output (MIMO) Node 

Now we consider the general case: a node with M input, N output links and C traffic commodities. Given 
are input demands S f and split ratios /3() specified per commodity, input priorities pi > 0, and output supply 
Rj. 


Define oriented demand: 

Sy = ftjSZ, 

and oriented priorities: 


Ef=i s\ 




Pij - Pi c 
Z^c= 1 


(3.15) 


(3.16) 


where i = 1,..., M, j = 1 ,... ,N and c = 1 ,... ,C. 

As before, we start by casting the allocation of input-output flows /,?• as a constrained optimization problem: 


M N C 

max 15] 41 > 


i—1 j =1 c—1 


subject to: 


(3.17) 


fij> 0, i = 1,..., M, j = 1,..., N, c = 1,..., C — non-negativity constraint; 
fij < Sfj, i = 1,...., M, j = 1,..., N, c = 1,..., C — demand constraint; 

M C 

4^ yy f?j < Rj, j = 1,..., N — supply constraint; 

*=1 C—l 

fii s 


E C fc' 

c' — l J ij Z^c' = l 


i — 1, • • •, M, j = 1,..., N, c = 1,..., C — proportionality 


constraint for commodity flows; 


E C xc 
c= 1 Ji' 


x^C £c x^N X^C c< 

2-jj=i 2^ c =i J ij L c =i 4 


C_1 , i = l, ...,M —FIFO constraint; 


(3.18) 

(3.19) 

(3.20) 


(3.21) 

(3.22) 


10 



(a) 


(b) 


For each input link i, such that 


spN v-^C fc 

Z-jj=l Z-/c=l J ij 


< Ef=i s h 


there exists an output link j 

Pi'j* Ef= 1 /«• > Pij* ES=i /f'j 


such that 
Vi' ^ i, 

and it is said that output j* restricts input flow i; 

For each input link i whose flow is restricted by some output j, 

Ef=i fu > 


E* 


Rj. 


priority constraint. (3.23) 


Here, the objective function (3.17) and constraints (3.18)-(3.21) are straight forward extensions of the objec¬ 
tive function (3.17) and constraints (3.2)-(3.5). The FIFO constraint (3.22) ensures that, for a given input 
link i, if any one of the output links is unable to accommodate its allocation of flow, all outflow from i is 
restricted proportionally (Daganzo, 1994, 1995). 


Let us focus the attention on the priority constraint (3.23). Condition (3.23)(a) says that if the flow from a 
given input link i is reduced by the output supply, there must be the most restrictive output link, which is 
denoted j*. If the output j* restricts flows from input links i' and i", then, according to condition (3.23)(a), 

C C 

Pi"J* 5Z fi'j* : Pi'3* 5Z fi"j* • 

c—1 c=l 

For each output link j, j = 1,... , N, flows Ylc= l fij f a b into two categories: (1) restricted by this output 
link; and (2) not restricted by this output link. Condition (3.23) (a) states that input-output flows of category 
1 are allocated proportionally to their priorities. For every output link j, j = 1,..., N, condition (3.23)(b) 
ensures that input-output flows of category 2 do not take up more output supply than was allocated for 
them when category 1 is non-empty. 


Recalling the MISO case, we can see that for iV = 1, formulation (3.17)-(3.23) translates to formulation (3.1)- 
(3.6). 


Remark. It must be noted here that when priorities pi are proportional to send functions X) c =i * = 
1,..., M, and M > 1, the condition (3.23) cannot be written simply as an extension of (3.7) for any output 
link j, j = 1,..., IV, as was done by Bliemer in (2007). He used the constraint 


E C £C 

C—1 j 1 j 

E C QC 
C= 1 °lj 


E C re 
c—1 J ij 

E C nc 

C— 1 °ij 


E C re 
c—1 J Mj 

EC QC 
C— 1 °Mj 


j = 1,..., N, 


(3.24) 


in the place of (3.23), which reduces the optimization problem (3.17)-(3.23) to an LR While suitable for 
N = 1, expression (3.24) is not an adequate representation of the constraint (3.23) for IV > 1, because it 
leads to the underutilization of the available supply, as will be shown in Section 3.3. 


Next, we present the algorithm for constructing the unique solution of the optimization problem (3.17)-(3.23). 


1. Initialize: 


o) 


Uj( o) 


k 


Rj 


0 ; 

i = 1,..., M, j = c = 1,... ,C. 


Uj(k) is the set of input links contributing traffic to output link j, whose input-output flows are still 
unassigned. 
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2. Define the set of output links that still need processing: 

If V(k ) = 0, stop. 

3. Check that at least one of the unprocessed input links has nonzero priority, otherwise, assign equal 
positive priorities to all the unprocessed input links: 


Pi(k ) = 


Pi, 


|Ujgv(fc) u i( k )\ ’ 


if there exists i' £ (Jjev(fc) ^j(k) ■ Pi' > 0, 
otherwise, 


(3.25) 


where 


Ujgy(fe) Uj( k ) denotes the number of elements in the union Uj<=V(fc) Uj(k)’, and for each output 


link j £ V(k) and input link i £ Uj(k) compute oriented priority: 

Pij{k) =Pi(k) ^ C c lS Z - 

2^ C =1 ^i 

4. For each j G V(k), compute factors: 

a(k) = 

A ’ e i&UjW Pij(ky 

and find the smallest of these factors: 

di*(k)= min aAk). 
jev(k) 

The link j* has the most restricted supply of all output links. 

5. Define the set of input links, whose demand does not exceed the allocated supply: 

r c \ 


(3.26) 


(3.27) 


(3.28) 


U(k) = lie Uj* (k) : £ s t < Pi(k)a r (k) \ . 


• If U ( k ) A 0i then for all output links j £ V ( k ) assign: 

ft, = Sf j , i£U(k), c=l 
__ c 

Rj(k + 1) = Rj{k)~ Y. E/b’ 
ieu(k) 0=1 

Uj(k + 1) = Uj(k)\U{k). 

• Else, for all input links i £ Uj»(k) and output links j £ V(k) assign: 

Pij(k)aj*{k) 


/ C _ QC 

ij 


E C q c f 

c'—l 


i £ Uj{k ) fl Uj• (A:); 


Rj(k+ 1 ) = Rj(k)- E E^; 

ieUj*(k) C=1 

Uj(k + 1) = Uj(k) \ Uj* (k), 
where, obviously, ?7j»(fc + l)=0. 

6. Set k := k + 1, and return to step 2. 


(3.29) 
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This algorithm takes no more than M iterations to complete. 

The following lemma states that in the case of N = 1, the MIMO algorithm produces exactly the same result 
as the MISO algorithm described in Section 3.1. 

Lemma 3.1. The MISO algorithm is a special case of the MIMO algorithm with N = 1. 


Proof. The proof follows from the fact that for N = 1, formulae (3.27)-(3.28) result in aj*(k) = aq(fc) = 


Ri(k) 

Etei/ijn) ft(fc) 


and j* = 1. 


□ 


The following theorem is the analog of Theorem 3.1. As in the MISO case, we leave the proof out, since a 
more general statement will be proved in Section 3.5. 

Theorem 3.2. The input-output flow computation algorithm constructs the unique solution of the maxi¬ 
mization problem (3.17)-(3.23). 


3.3 Relationship with Tampere Et Al. and Bliemer Node Models 

It can be shown that the node model described here is a slight generalization of the Tampere et al. unsignal¬ 
ized node model developed in Tampere et al. (2011), In particular, for a specific choice of priorities, Pi, the 
equations of section 3.2 become exactly the equations of this model. 

Theorem 3.3. By setting the priorities for the MIMO model equal to the link capacity, F, 

Pi — -E, 

with Fi the capacity of link i, we obtain the Tampere et al. node model. 


Proof. The Tampere et al. model does not deal with multiple commodities, so first, let C = 1 and drop the 
commodity index c. Then, our oriented priorities (3.26) become: 


Pij = Ft 


Si ' 


Again for simplicity, assume that we do not have any input links whose demands can be wholly met. Then, 
the smallest Oj(fc) factor is 


Ri 


mm - 


(k) 


Fi 


Si 


and the resulting flows are 


fij = Si 


Pij{k)aj* (k) 
' o.. 


= p ij (k)a J 


■W = «|i 


Ri 


E l 7 1 

iGUi* S, 


(3.30) 


comparing the above development to equations (19), (24), (26), and (28) in Tampere et al. (2011), we have 
recovered exactly the Tampere et al. node model. □ 
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Figure 2: A 4-by-4 example intersection. 


Example. This specific example is taken from Tampere et al. (2011), to illustrate both the procedure of our 
MIMO node algorithm and the equivalence of our node model to the Tampere et al. model proved above. 

Consider the 4-by-4 intersection shown in Figure 2. Suppose that our relevant parameters are: 


C= 1 
Si = 500 
Pi = F\ = 1000 
R 5 = 1000 
# 15 = 0 
#25 = 0.05 
#35 = 0.125 
#45 = V 17 


S 2 = 2000 
p 2 = F 2 = 2000 
R 6 = 2000 
#16 = 0.1 
#26 = 0 
#36 = 0.125 

#46 = 8 A 7 


S 3 = 800 

p 3 = f 3 = 1000 

R 7 = 1000 
#17 = 0.3 
#27 = 0.15 

#37 = 0 
#47 = 8 / 17 


S 4 = 1700 

Pa = Fa = 2000 
R s = 2000 
#18 = 0.6 
#28 = 0.8 
#38 = 0.75 

#48 = 0 


Our solution algorithm proceeds as follows: 


First, compute oriented demands: 


Sis — 0 

S 25 = 100 
S 35 = 100 
S 45 = 100 

k = 0 : 


Si 6 = 50 
S26 = 0 

s 36 = 100 
S 46 = 800 


S n = 150 
S 27 = 300 

S37 = 0 
S 47 = 800 


Si 8 = 300 
S 2s = 1600 
S 38 = 600 
S48 = 0. 


1 . U 6 (0) = {2,3,4}; C/ 6 (0) = {1,3,4}; U 7 (0 ) = {1,2,4}; U 8 ( 0) = {1,2,3} 

2. E(0) = {5,6,7,8} 

3. No adjustments to the pi are made. The oriented priorities are: 

#15 = 0 #16 = 100 #17 = 300 #18 = 600 

#25 = 100 #26 = 0 #27 = 300 #28 = 1600 

#35 = 125 #36 = 125 #37 = 0 #38 = 750 

#45 = 118 #46 = 941 #47 = 941 #48 = 0 


4. a 5 (0) = 2.92; a 6 (0) = 1.72; a 7 (0) = 0.649; a 8 (0) = 0.678 

a r ( 0) = a 7 (0) = 0.649 

5. 17(0) = {1}, as Si = 500 < #i(0)a 7 (0) = 1000 x 0.649 

• fi6 — Si 6 = 50; /i 7 = Si 7 = 150; /is = Sis = 300; 

• R 6 ( 1) = 2000 - 50 = 1950; R 7 ( 1) = 1000 - 150 = 850; R 8 ( 1) = 2000 - 300 = 1700 
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• u 5 ( 1 ) = {2,3,4}; U$(l) = {3,4}; U 7 ( 1 ) = {2,4}; U 8 ( 1 ) = {2,3} 

k = 1 : 

2. 17(1) = {5,6,7,8} 

3. No changes are made to priorities or oriented priorities. 

4. a 5 (l) = 2.92; a 6 (l) = 1.83; a 7 (l) = 0.685; a 8 (0) = 0.723 

aj.(l) = a 7 (l) = 0.685 

5. 17(1) = 0, as S 2 = 2000 ^ p 2 (l)a 7 (l) = 1370, and S 4 = 1700 ^ p 4 (l)a 7 (l) = 1370, 

• / 25 =P2b(1)ot(1) = 68.5 

• hr = p2 T (l)ar(l) = 205.5 

• hs = P 2 8 (l)ar(l) = 1096 

• / 45 — ^ 45 ( 1 ) 0 . 7 ( 1 ) = 80.6 

• he = p 4 6(l)a r (l) = 644.5 

• / 47 = p 47 (l)a r (l) = 644.5 

• £ 5 ( 2 ) = 1000 - 68.5 - 80.6 = 850.9; R 6 { 2) = 1950 - 644 = 1305.5; 

R 7 ( 2) = 850 - 644.5 - 204.5 = 0; # 8 ( 2 ) = 1700 - 1096 = 604 

• C/ 5 (2) = {3}; U 6 ( 2) = {3}; !7 7 (2) = 0; J7 8 (2) = {3} 

k = 2 : 

2. 17(2) = {5,6,8} 

3. No changes are made to priorities or oriented priorities. 

4. 015 ( 2 ) = 6.81; a 6 (2) = 10.45; a 8 (2) = 0.805 

a r { 2) = a 8 (2) = 0.805 

5. U( 2) = {3}, as 5 3 = 800 < p 3 ( 2 )a 8 ( 2 ) = 805 

• /35 = S'ss = 100; fs6 — S 36 = 100; / 38 = S 38 = 600; 

• Uj( 3) = 0 Vj 


fc = 3 : 

2. 17(3) = 0, so the algorithm terminates. All flows have been found. 

As discussed in Section 3.2, for N > 1 the demand-proportional priority node model of Bliemer (2007) 
cannot be written as a special case of this node model, and that replacing the priority constraint (3.23) with 
a constraint of the form (3.24) leads to underused supply. Let us consider a simple example: 

Example. Suppose we have a 2-by-2 intersection, with parameters 

<7=1 pi = 0.5 p 2 = 0.5 

Si = 1000 S 2 = 1000 Ri = 600 R 2 = 1000 

Pll = 9 / 10 /?12 = V 10 

/?21 = 0 /? 22 = 1 
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It should be clear that our algorithm will solve this problem in two iterations. In the first iteration, we will 
have 


Ol (0) = 600/goo = 2 /3 = Oj.(0) 

02 ( 0 ) = 1000 /1100 = 10 /ll 

and, as U 3 * (0) = t/i(0) = {1}, the flows found will be 
fu = 600 

hi = 2/3 x 100 « 67. 

In the second iteration, link i = 2 will fill up all remaining space in link j = 2: 

hi = 0 

hi =Ri~ hi » 933, 

and the total flow passing through the node is 
fu + hi + hi + hi = 1600. 


On the other hand, in the Bliemer model we first encounter the supply restriction Ri that, together with 
the FIFO constraint (3.22), reduces the flow from input link 1 by the factor 2 / 3 , and we get 

hi = 600 

fu = 2 /3 X 100 « 67, 

just as in our calculation above. There is no flow from input link 2 to output link 1: hi = 0. Next, 
condition (3.24) implies that 

fu = hi + hi = 600 
S 11 Sia + S 22 900 

which yields: 

hi = ^(<Si2 + Sii) ~ hi ~ 667, 
and so, the total flow through the node is 
/11 + /12 + /21 + fii = 1334. 

As we see, the Bliemer model violates requirement 2 from the list in the beginning of Section 3: it does 
not maximize the total flow through the node. Failure to maximize the flow happened in the LP setting 
of Bliemer (2007) due to the incorrect constraint (3.24) that was imposed in pursuit of an easy formulation 
of the optimization problem. 

It was shown in Tampere et al. (2011) that the Bliemer model also violates the invariance principle — 
requirement 7 from the list. The violation of the invariance principle happened because supply allocation 
was governed by the input demand. The input link started to congest, and the demand jumped in infinitely 
small time, which immediately triggered change in the supply allocation of the output link, ultimately 
resulting in discontinuous flow. In a discrete time system, such as LNCTM, these considerations may be 
simply ignored. However, since our proposed node model is generic and can be used outside of the LNCTM, 
we should note that to avoid violating the invariance principle, input link priorities must not depend on 
traffic state. 
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3.4 Single-Input-Multiple-Output (SIMO) Node with Relaxed FIFO Condition 


Sometimes the FIFO rule may be too restrictive. An example of such situation is a junction, where a 1-lane 
off-ramp diverges from a 5-lane freeway, especially, if this off-ramp is relatively short and terminates with 
a signal. Jammed off-ramp blocks the flow on the whole freeway, which is not realistic. A more realistic 
behavior would be for this jammed off-ramp to only partially restrict the mainline flow. For instance, one 
could say that this off-ramp restricts only 40% of the mainline traffic. On the other hand, if the mainline 
happens to be jammed, we want 100% of the off-ramp traffic to be restricted. Thus, in a given diverge 
junction FIFO restrictions on outgoing flows are not symmetric. The model of a diverge node with a partial 
FIFO restriction is described next. 


Consider a node with a single input and N output links. Given are the input demand per commodity, 
S c , c = 1 ,C, the output supply Rj, split ratios •, j = 1,..., N, and the mutual restriction matrix 
£ Hl NxN , r/jji £ [0,1], ?jjj = 1. Element rjjj> of the mutual restriction matrix specifies the portion of 
the flow in the output link j' affected by the restriction of the output link j. In the example above, where 

the mainline output link is identified as 1 and the off-ramp as 2, this matrix assumes the form ^ j 

To formulate the optimization problem for the SIMO case, we can re-use the objective function (3.17) and 
constraints (3.18)-(3.21) directly, with M = 1; constraint (3.23) can be dropped, since there is no competition 
for the output supply between incoming flows; and the FIFO constraint (3.22) has to be replaced. 


Using (3.15), we obtain oriented demand: 

i = 1,... ,7V. 


Q c - 
Nj - , 


The relaxed FIFO constraint can now be written 
c c 

Y fij < a - vj'j) Y s h+ 



, j = !,•••, N, V/ ^ j. 


(3.31) 


(3.32) 


The right hand side of the last inequality consists of two terms. The first term represents the portion 
of the oriented demand unaffected by a possible supply shortage in the output link j'. If the full FIFO 
rule is enforced, that is, rjj'j = 1, then this term equals 0. The second term represents the portion of the 
oriented demand influenced by the output j'. The necessary, but not sufficient, condition for activation of 


this constraint is 
term equals 0. 


I2c=i si 


< 1 for some j'. If the FIFO rule is abandoned, that is, rjj’j = 0, then the second 


The algorithm for computing input-output flows /f • follows. 


1. For each output j, compute the reduction factor: 


Rj 

a.j = mm < 1. — 7 = r ^ - 

l Eti sij 


(3.33) 


The factor ctj defines by how much the demand directed to the output link j must be scaled down to 
satisfy the supply constraint Rj. 


2. For each output j, compute the total flow directed to this output link: 
hi = min j otj Y S ^v fgj \ “ 'Vj'j) Y S h + Vj'jaj' Y S h \ 

l C=1 l C—1 c—1 J . 


(3.34) 


Here we pick the term that restricts the flow to output link j the most: the restriction may come from 
the output link j itself, or from other output links, denoted j', that influence link j through mutual 
restriction coefficients rjj'j. 
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(3.35) 


3. Assign input-output flows per commodity: 

f» = Sii J' 1 , j = 1-JV. 

2^ c '=i ‘-’lj 

Example. Suppose we have one input and three outputs, with parameters: 

(7=1 S = 1000 
i?i = 100 R 2 = 400 R 3 = 300 

fti = 0.2 ft 2 = 0.5 ft 3 = 0.3 

V 12 = 0.2 7721 = 1 7731 = 0 

V13 = 0 7723 = 0.5 7732 = 0 

Our solution algorithm for finding the resulting flows proceeds as follows: 

1. fti = 200; ft 2 = 500; ft 3 = 300; 

„ 100 „ r 400 300 , 

2 - ai = 200 = °- 5; “ 2= 500 = °- 8i as = 300 = 


3. 


/n = min 0.5 • 200, min (1 - 1)200 + 1 • 0.8 • 200, (1 - 0)200 + 0-1-200 


= min i 100, min < 160, 200 

fn = 100 

/12 = min i 0.8 • 500, min \ (1 - 0.2)500 + 0.2 • 0.5 • 500, (1 - 0)500 + 0-1-500 


= min ^ 400, min <, 450, 500 

fi2 = 400 

/i 3 = min <| 1 • 300, min { (1 - 0)300 + 0 • 0.2 • 300, (1 - 0.5)300 + 0.5 • 0.8 • 300 


= min i 300, min < 300, 270 

/is = 270 


Finally, we state the result of this Section as a theorem. 

Theorem 3.4. The SIMO input-output flow computation algorithm constructs the unique solution of the 
maximization problem (3.17)-(3.21), (3.32) with M = 1. 


3.5 MIMO Node with Relaxed FIFO Condition 


For a node with M input and N output links we are given input demands Sf and split ratios ftj specified per 
commodity, input priorities pi > 0, output supply Rj and mutual restriction matrices {r)jj'} % (v’jf G [0,1], 
rfjj = 1), i = 1,..., M, j = 1 ,,N and c = 1,..., (7. Oriented demand S+ and oriented priorities p l: j are 
defined in (3.15). 


To formulate the optimization problem for the MIMO case with the relaxed FIFO constraint, we extend the 
SIMO relaxed FIFO constraint (3.32) to arbitrary M: 




C— 1 


C— 1 


( EtifrA 

VEf=i S9.,J 



, i = j = (3.36) 
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If = 1 for all i, j and j', then constraint (3.36) is equivalent to (3.22). 

The optimization problem (3.17)-(3.21), (3.23), (3.36) satisfies the requirements listed in the beginning of 
Section 3. The following algorithm for computing input-output flows /() generalizes the MIMO algorithm 
from Section 3.2 taking into account constraint (3.36). 


1. Initialize: 


^( 0 ) 

:= R: 

^•(0) 

■■= u j; 

sm 

- CC . 

’ 

C 

4(o) 

■■= E4(°); 

C— 1 

k 

:= 0 ; 


2. Define the set of output links that still need processing: 


V(k) = {j : Uj(k) ± 0} . 

If V{k) = 0, stop. 

3. For each input link i £ Ujey(fe) Uj(k), calculate input link priority Pi(k) according to expression (3.25). 

4. For each output link j £ V(k) and input links i £ Uj(k), calculate oriented priorities: 

Pij(k) = Pi{k)— — (3.37) 

2^j'ev(k) bij'{k) 


5. For each j £ V(k), compute factors a,j(k) according to (3.27) and find the most restrictive output link 

6 . Define the set of input links, whose demand does not exceed the allocated supply: 


U(k) = \ i £ Uj * (k) : 


Sii( k ) ) <Pi(k)aj*(k) 

yjev(k) 


If U ( k ) ^ 0, then for all output links j £ V ( k ) assign: 

ftj = SUV, i€U(k), c = 1 ,..., C; 


Rj (fc + 1) 

Uj(k + 1) 


E 

ieir(k) c=1 


Uj{k) \ U(k). 


(3.38) 
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• Else, for all input links i £ Uj * (k), output links j £ V ( k ) and commodities c = 1,..., C, assign: 


%•(* + !) 

where 


Sij(k + 1 ) 


Sfj{k + 1) 
Rj(k + 1) 


sm 


Pij{k)aj- (k) 

E £=!%•(*)’ 


<%•(*;) 


Sij{k + 1 ) 


J 


j : rf rj = 1; 

Vj*j < 1, « e Uj(k) n Uj*(k); 


min ( Sij(k), (1 - rfj.j) ^ S c zj + 1 ^ 

2 ^c=i 


C=1 


Su(k), itu J (k)nu j .(ky, 


c 


U k )~ E E^; 

ieUj* =1 c=l 



(3.39) 

(3.40) 


(3.41) 


Uj(k + 1) = Uj(k) \{i£ Uj- (k) : r\),j = l} . 

Here, the assignment (3.39) reproduces that of the original MIMO algorithm, (3.29), and the 
expression (3.41) represents the selection of the minimum in the SIMO algorithm, (3.34), where 
mutual restriction coefficients {T)j*j}‘ are applied. Equation (3.40) is an application of our pro¬ 
portionality constraint for commodity flows, (3.21). 


7. Set k := k + 1, and return to step 2. 


This algorithm takes no more than (M + N — 2) iterations to complete. 

Remark. An intuition for the composition in (3.41) may be expressed as, “is the currently-considered output 
j* the most restrictive output on the movement (i. j) (RHS of the minimum), or has j already had a greater 
restriction enforced upon it (LHS of the minimum)?” 

The following lemma states that in the case M = 1, the MIMO algorithm with relaxed FIFO condition 
produces exactly the same result as the SIMO algorithm with relaxed FIFO condition described in Section 3.4. 

Lemma 3.2. The SIMO algorithm with relaxed FIFO condition is a special case of the MIMO algorithm 
with relaxed FIFO condition with M = 1. 


Proof. In the case of M = 1, factors aj(k), defined in (3.27), reduce to: 


;; ,m 'bEE-SU* 

3 Oc ( U\ 

E c =i s ij( k > 


Thus, 


Rj Ec=i s i( k ) 

j = arg min = arg min 


i?, 


£Ei SUk) 


Ec=i W 


Ri* 


E“ =1 sj,,(fe) 


< 1. 


MIMO algorithm with M = 1 goes into iteration k+1 only if U ( k ) = 0, which is equivalent to 

Then the assignment of input-output flows /f • for output links j, such that ryl. • = 1, is given by (3.39) that 
translates to: 


f'f. = a j -{k)SlAk), c = 1,..., C, 


where 


ij-(k) = 


Rj * 


E C qc 

C— 1 °lj* 


= min < 1, 


Rj* 


Ec=l S lj* ( k ) 
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which is the same as (3.33) for the output link j*. Obviously, Xl c =i fij = fij satisfies (3.34) for all j, such 
that = 1, which always includes j*. For output links j, such that ijj*j < 1, formula (3.41) achieves the 
same result as (3.34) through the iterative pairwise comparison. □ 


The next lemma states that the when FIFO condition is fully enforced, that is, rf^y = 1 for all output link 
pairs {j. j'} and all input links i, the MIMO algorithm with relaxed FIFO condition reduces to the original 
MIMO algorithm from Section 3.2. 

Lemma 3.3. The original MIMO-with-FIFO algorithm is a special case of the MIMO algorithm with relaxed 
FIFO condition with mutual restriction coefficients r]F, = 1, i = 1,..., M, j,f = 1,... ,N. 

Proof. The proof follows from the fact that steps 6 of the MIMO-with-FIFO and the MIMO-relaxed-FIFO 
algorithms are equivalent when rffy = 1 (compare expressions (3.29) and (3.39)). □ 

Remark. Note that in the case of multiple input links mutual restriction coefficients are specified per input 
link. This can be justified by the following example. In the node representing an intersection with 2 input 
and 3 output links, shown in Figure 3, consider the influence of the output link 5 on the output link 4. If 
vehicles enter links 4 and 5 from link 1, then it is reasonable to assume that once link 5 is jammed and 
cannot accept any vehicles, there is no flow from 1 to 4 either. In other words, r]l 4 = 1. On the other hand, 
if vehicles arrive from link 2, blockage of the output link 5 may hinder, but not necessarily prevent traffic 
from flowing into the output link 4, and so ?y| 4 < 1. 


A 

5 


4 

< 

1 


A 


2 


3 




Figure 3: Intersection node. 

Remark. As a rule of thumb, we suggest setting the restriction coefficient rfy equal to the portion of lanes 
serving movement that also serve movement As an example, if movement {i,f) is served by four 

lanes, and of those lanes, one serves movement (i,j), then a reasonable value for is 0.25, as blockage of 
movement (i. j) will block 25% of (f,j')’s lanes. 

Example. Recall our earlier example for a 4-by-4 junction during our discussion of the MIMO-with-FIFO 
node model in Section 3.2. Recall that our example supposed that the input links (see Figure 2 for the 
network topology) had capacity values of F\ = F 3 = 1000, and F 2 = F 3 = 2000. Say that the links with 
capacity of 1000 are one-lane roads, while those links with capacity of 2000 are two-lane roads, with left and 
right turns permitted from one lane each, and the through movement permitted from both lanes. Following 
our rule of thumb, this would lead to a set of reduced mutual restriction coefficients: 

? ?76 = ? ?56 = 0-5 V 75 = V 57 = 0 
? ?58 = ? ?78 = 0-5 V 75 = ^57 = ^ 

with all other restriction coefficients taking a value of 1. Let us examine the effect of our FIFO relaxation on 
this example. Since C = 1, we ignore the index c and perform some cancellations to clean up the notation. 

k = 0 : This iteration will proceed exactly as before, with U{ 0) = {1} and flows from link 1 taking the 
demand value. 

k= 1 : 
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2 . V(l) = {5,6,7,8} 

3. No changes are made to priorities for any i £ Ujevfi) 

4. No changes are made to oriented priorities. 

5. a 5 (l) = 2.92; a 6 (l) = 1.83; a 7 (l) = 0.685; a 8 (0) = 0.723 

aj.(l) = a 7 (l) = 0.685 

6 . 17(1) = 0, as S 2J {1) = 2000 £ pa(l)a 7 (l) = 1370 
and 5 4j -(1) = 1700 ^ p 4 (l)a 7 (l) = 1370 

• Ur = p 47 (l)a 7 (l) = 644.5 

• Hr = P 2 7 (l)a 7 (l) = 205.5 

• $28(2) = (1 ~ Vrs)S28 + ?? 7 g x ^ 27 /S27 x $28 

= (1 - 0.5) x 1600 + 0.5 x 205.5/300 x 1600 
= 1348 

• # 46 ( 2 ) = (1 - 7 ? 76 )$ 46 + ?? 76 X /47/S47 x $46 

= (1 - 0.5) x 800 + 0.5644.5/goo x 800 
= 722.25 

• R 7 ( 2) = 850 - 644.5 - 204.5 = 0 

• J7 e (2) = {2,3,4}| C/ 6 (2) = {3,4}; J7 r (2) = 0; C/g(2) = {2,3} 


k = 2 : 


2. V(2) = {5,6,8} 

3. No changes are made to priorities pi. 

4. Due to changes in two oriented demands, the new oriented priorities for input links 2 and 4 must be 
computed: 

• p 25 (2) = 20 00 x 100/1448 = 138.12 

• p 28 (2) = 2000 x 1348/4448 = 1861.9 

• P 45 (2) = 20 00 X 100/822.25 = 243.24 

• p 46 (2) = 2000 X 722.25/822.25 = 1756.8 

5. The coefficients a ; (2) are calculated with the new oriented priorities: 

. „ (n\ _ _ fis( 2 ) _ _ _ 1000 _ _ 1 09 

“S'- ' P25(2)+P35(2)+P4 5 (2) 138.12+125+243.24 

. „ / 9 \ _ ^6 (2) _ 1950 _ i n1 

0,6 +P36(2)+p 46 (2) ~ 125+1756.8 “ i ' U ‘ i 

• as(2) = +p 2S (2)+^ 38 (2) = i86i 7 9°+750 = 0-651 = %*( 2) 

6. 17(2) = 0, as S 2j { 2) = 1448 £ p 2 (2)a 8 (2) = 1301 
and Ej $3j(2) = 800 £ p 3 (2)a 8 (2) = 651, 

• Hs = p 2 8(2)a 8 (2) = 1861.9 X 0.651 = 1211.75 

• fas = p 3 8(2)a 8 (2) = 750 X 0.651 = 488.25 

• nl b = 1, so / 25 = p 2 5(2)a 8 (2) = 138.12 X 0.651 = 89.916 

• r?| 5 = 1, so / SB = p 35 (2)a 8 (2) = 125 X 0.651 = 81.375 
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• ryle = 1, so / 36 = p 36 (2)a 8 (2) = 125 X 0.651 = 81.375 

• R 5 ( 3) = 1000 - 89.916 - 81.375 = 828.71; R 6 ( 3) = 1950 - 81.375 = 1868.6 

Rsl 3) = 1700 - 1211.75 - 488.25 = 0 

• U B ( 3) = {4}; Uq(3) = {4}; U 8 ( 3) = 0 


2. U(3) = {5,6} 

3. No changes are made to the priorities p 4 , those of the only link i remaining. 

4. The oriented demands of link 4 remain the same as calculated in the previous iteration. 

5. a 5 (3) = 828.71/258.39 = 3.21 

a 6 (3) = 1861 . 6 / 1741 . 6 = 1.073 = a r (3) 

6. U( 3) = {4}, as Sij = 774 < p 4 (3)« 6 (3) = 2146 

/ 45 = ^45(3) = 100; / 40 = S 40 (3) = 722.25 

Uj( 4) = 0 Vj 


2. U(4) = 0, so the algorithm terminates. All flows have been found. 


Comparison of the resulting flows with those from the non-relaxed FIFO example illustrates several differences 
caused by this relaxation. Unsurprisingly, the two-lane links (2 and 4) that take advantage of the FIFO 
relaxation send greater amounts of vehicles to links 5, 6 and 8 . Previously, these flows had been cut off by 
the filling of link 7’s supply. Perhaps more interesting is the effect on the flows of link 3. In the previous 
example, link 3 was able to send vehicles equal to its demand, since there was enough leftover supply in link 
8 after the flow / 2 s was constrained by FIFO. In our example with relaxed FIFO, however, link 8 had more 
supply consumed by the higher-priority link 2, and it was found that link 3 was unable to send its entire 
demand. Since we had assumed that link 3 was a one-lane road and thus bounded by FIFO, this restricted 
its other flows as well. 

In both cases, link 3 was lower priority than link 2 and was disadvantaged in claiming link 8 ’s supply, but 
qualitatively, permitting a FIFO relaxation for the two-lane roads exacerbated the “advantage” they already 
had over the one-lane roads (from having priority equal to capacity). 

The main result of this Section can be stated as the following theorem. 

Theorem 3.5. Given a set of input links i £ 1,..., M, output links j 6 1..Y, commodities c € 1,..., C, 

priorities pt, split ratios /?,, and mutual restriction coefficients { 77 *• •,}, the algorithm of Section 3.5 obtains 
the unique solution of the optimization problem (3.17)-(3.21), (3.23), (3.36). 


Proof. The priority constraint (3.23) makes this optimization problem non-convex, except in the above- 
mentioned special case where N = 1 and the priorities pi are taken as proportional to the send functions. 
Recall our discussion in Section 3.2 that the priority constraint must be replaced entirely with the unrealistic 
constraint 3.24 to obtain an LP. Our priority constraint increases the complexity of the problem significantly. 
In fact, we conjecture that in the MIMO case, both with and without relaxation of the FIFO constraint, 
there is no way to verify a solution faster than re-solving the problem (3.17)-(3.21), (3.23), (3.36). Thus, we 
will prove optimality by showing that as our algorithm proceeds through iterations, it constructs the unique 
optimal solution. 


23 



We may decompose the problem into finding the M ■ N ■ C interrelated quantities {/£■}. The C flows for 
each (i,j) are further constrained by our commodity flow proportionality constraint (3.21); solving for one of 
{flj, ■ ■ fi'j} also constrains them all. Our task then becomes finding optimal values for each oiM-N subsets 
{ flj, ■ ■ ■, fij}- Our algorithm finds at least one of these M ■ N subsets per iteration k. The assignments 
are done by either equation (3.38) or equation (3.39). Over iterations, subsets are assigned to build up the 
unique optimum solution. We can show that each subset assigned is optimal; that is, at least one of the 
constraints is tight. 


Let us first consider formulae (3.40)-(3.41), our implementation of the relaxed FIFO constraint. In step 5 of 
our algorithm, we identify a single output link as j*. The argmin in this step picks out a single link as the 
most restrictive of all output links. By the relaxed FIFO construction, all i £ Uj* will feel a relaxed FIFO 
effect instigated by j* as the most restrictive link. For a generic j A j*, equation (3.41) enforces the relaxed 
FIFO constraint by decaying the oriented demand S t j. In fact, this oriented demand acts as a proxy for 
the restricted FIFO constraint. Since a /,?■ that is restricted by relaxed FIFO will, by construction, never 
obtain /() = .S')), the running quantity Sfj(k) represents an “effective” oriented demand after application of 
restricted FIFO. If it turns out that a flow is found ff- = SfAk ) < for some k, then this flow value has 
an active constraint imposed by relaxed FIFO. 


Now consider flows assigned by equation (3.39). In the special case where j = j *, we have 

t E /«•=£ E 4 -«=)|sM- 

c=l ieu r (fc) C=1 ieu^(k) a ij* \ K ) 

= Pij-(.k)a r (k) 

ieu^(k) 


= Pij-(k) 

yieu^ik) 

= Rj*(k), 


RAk) 


Ei 


ieUj*{k) 


Pij * ( k ) 


so all the flows into link j* take up all available supply. These flows are thus constrained by the supply 
constraint, (3.20). 

For flows assigned by equation (3.39) in the case of j A 1 J*> w e have i £ Uj*, and a restriction coefficient 
rjj*j = 1 that implies a strict FIFO constraint. In other words, we need 

EtiA __ EtiW 

i (fij + fij•) Eti (% (k) + S£, (A;)) 

Examining the denominator of the LHS, we can write 


C 


A + A * = ( Pij ( k ) +Pij*{k)) cij*(k) 


C— 1 

C 

= £■ 


Pi( k ) 


Ejev(k) S ij( k ) 
Straightforward algebra then gives 


(Sij(k) + 5y*(fc)) a,j*{k). 


E C £C 

c—1 


£c=i (fij + ffj *) 


Pi(k) ^ Sjj(k) ^ J2jev(k) Rjj(k) ^ gj* (fc) 

A[Pi{k) Sfj (k) + Sjj, (k) Y,jev { k)Sm ^.(k) 
Eti Sfj(k) 

Eti (S?j(k) + s?A k ))' 
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as originally required. Thus, flows {//)}, j ^ j* assigned in equation (3.39) are constrained by the strict 
FIFO constraint, (3.36) with = 1. 

For flows assigned by equation (3.38), we have for each flow i, f9. = SA(fc), so these flows are constrained 
by (3.19). 


Uniqueness of the solution is equivalent to showing that we must assign flows in the order described. That 
is, that flows for the most restrictive link j* should be assigned at iteration k. It turns out that setting the 
flows of any other j ^ j* may produce solutions that violate the problem constraints. At iteration k, let j' 
be some other constrained outgoing link, with aj * (fc) < ay (fc), and let i' be any input link in Uj * (k) (~l Uj> (k). 
Proceeding through our algorithm, we will obtain an upper bound on the sum of flows from i' to j' through 
the restriction coefficient: 

G 

fi'j < < min | Si'j>(k + 1 ),pi<{k + l)ay{k + 1) j 

C— 1 


— Si’j 1 {k + 1) 


Vi* 


C— 1 


E C £C 

c= 1 J i'j* 




Y. s h = k. 

C—1 


Suppose instead that we had decided to set flows into j' at iteration k. Then we would have 

C 


yfh =Pi'j'(k)a j »(k) = ^ n ^ Rji(k). 


C— 1 


E 


ieu^ik) Pi'j' 


Vi'o' (k)' 


(3.42) 


Note that (3.42) does not make use of the restriction coefficient and therefore does not take (relaxed) 

J2 C - f c i * 

FIFO into consideration. In fact, for a suitably large rj), •, and/or small ratio , we can obtain 

' ' ' b c =!*,/,•» 


Pi’j'jk ) 

Eieu^i^Pi'j'ik) 


«»■(*)> k Ey 1 ' +i-4r 

\ 2-^c= i ^i'j* / 


C 

E'% 


J V3 


fi'i> > fv. 


c—1 

breaking our upper bound imposed by (relaxed) FIFO (3.36). 

Corollary 3.1. Theorems 3.1, 3.2, and 3-4 follow from the Theorem 3.5 as special cases. 


□ 


Remark. It was discussed above that the Tampere et al. node models are special cases of this node model 
(with all ijfy = 1 and certain choices of pf). The above proof then also applies to the similar theorem for 
that model. 


4 Node with Undefined Split Ratios: Traffic Assignment 


In this section we consider a MIMO node with M input links, N output links and C commodities, where 
some of the split ratios /?,■) are not defined a priori and must be computed as functions of the input demand 
Sf, priorities pi and the output supply Rj, i = 1,..., M, j = 1,..., N and c = 1,..., C. This may occur 
if the modeler chooses to let drivers at a node select a route to their destination in response to changing 
conditions; a specific example of such a node is given in an example describing an interchange between 
a freeway and a managed lane to be discussed after presentation of the algorithm. Here we present the 
algorithm for computing undefined split ratios based on the following definitions and assumptions: 
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• Define the set of commodity movements, for which split ratios are known, B = c /%€[ 0,1]}, 

and the set of commodity movements, for which split ratios are to be computed, B = 
{{*, j, c} : 0% are unknown}. 

• For a given input link i and commodity c such that Sf = 0, assume that all split ratios are known: 
{i,j,c} £ B. 9 

• Define the set of output links, for which there exist unknown split ratios, V = (j : 3{i,j, c} £ 6}. 

• Assuming that for a given input link i and commodity c split ratios must sum up to 1, define the 

unassigned portion of flow = 1 - ffij- 

• For a given input link i and commodity c such that there exist {i,j,c} £ B , assume f3 i > 0, otherwise 
the undefined split ratios can be trivially set to 0. 


• For every output link j £ V, define the set of input links that have an unassigned demand portion 
directed toward this output link, Uj = {* : 3 {i. j, c} £ B }. 

• For a given input link i and commodity c define the set of output links, split ratios for which are to be 
computed, Vf = {j : 3i £ Uj}, and assume that if nonempty, this set contains at least two elements, 
otherwise a single split ratio can be trivially set to /3}. 

• Assume that input link priorities are nonnegative, Pi > 0, i = 1,... ,M, and Pi = 1- Note that, 
although in section 3 we did not require the input priorities to sum to one, in this section we assume 
this normalization is done to simplify the notation. 


• Define the set of input links with zero priority: U zp = {i : Pi = 0}. To avoid dealing with zero input 
priorities, perform regularization: 


Pi = Pi 



j l_\UzA = M-\U zp \ 
+ M M Pi M 


Wzp\ 
M 2 ’ 


(4.1) 


where \U zp \ denotes the number of elements in set U zp . Expression (4.1) implies that the regularized 
input priority pi consists of two parts: (1) the original input priority pi normalized to the portion of 
input links with nonzero priorities; and (2) uniform distribution among M input links, jj, normalized 
to the portion of input links with zero priorities. 

Note that P; > 0, * = 1,..., M, and i Pi = ' ■ 


The algorithm for distributing j3\ among the commodity movements in B , that is assigning values to the 
a priori unknown split ratios, aims at maintaining output links as uniformly occupied as possible. It is 
described next. 


1. Initialize: 


W) 

Tm 

um 

no) 

k 


fa 

Uj ; 

13 ; 

o, 


if £ B, 

otherwise; 


Here Uj(k) is the remaining set of input links with some unassigned demand, which may be directed 
to output link j; and V(k) is the remaining set of output links, to which the still unassigned demand 
may be directed. 

®If split ratios were undefined in this case, they could be assigned arbitrarily. 
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2. If V(k) = 0, stop. The sought-for split ratios are j/3,?(/c) j, i — 1,..., M, j = 1,... ,N, c = 1,... ,C. 

3. Calculate the remaining unallocated demand: 

3 c i (k) = F i (k)Sf, i = c = 1,... ,C. 


4. For all input-output link pairs calculate oriented demand analogously to (3.15): 

S?j(k) = feWSf. 


5. For all input-output link pairs calculate oriented priorities: 


Pij{k) 

- ES=i7 ZjSf 

= Pi 9C 

Z^ C =1 


(4.2) 

with 




7 ij(k) 

j ftp 

\ 4( fc ) + m’ 

if split ratio is defined a priori: {i,j, c} G B , 
otherwise, 

(4.3) 


where |V) C | denotes the number of elements in the set V['. Comparing the expression (4.2)-(4.3) 
with (3.26), one can see that split ratios /JG(fc), which are not fully defined yet, are complemented 
with a fraction of /3^(k) inversely proportional to the number of output links among which the flow of 
commodity c from input link i can be distributed. 

Note that in this step we user regularized priorities pi as opposed to the original pi, i = 1,..., M. This 
is done to ensure that inputs with Pi = 0 are not ignored in the split ratio assignment. 


6. Find the largest oriented demand-supply ratio: 


p+(k) 


= max max 
j * 


Pij ( k)Rj 


E p^ k )- 

ieUj 


7. Define the set of all output links, where the minimum of the oriented demand-supply ratio is achieved: 


Y(k) = arg min min 

jev(k)ieUj{.k) 


J2c=i Sfjik) 

Pij(k)Rj 


E 

isl¬ 


and from this set pick the output link j with the smallest output demand-supply ratio is minimal 
(when there are multiple minimizing output links, any of the minimizing output links may be chosen 
as j 


7 = arg mm 

jer(k) 


s^k) 

R; 


8. Define the set of all input links, where the minimum of the oriented demand-supply ratio for the output 
link j~ is achieved: 


Wj- ( k) = arg min 


Ec=i stUk) 


ieu,-(k) Pij- (k)Rj- 


E 

ieu.- 


Pij-(k), 


and from this set pick the input link i and commodity c with the smallest remaining unallocated 
demand: 

{i~,c~} = arg min S°(fc). 
i G Wj-(k), 
c : X- (*:) > 0 
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9. Define the smallest oriented demand-supply ratio: 


P (fc) = 


E c =i S?...(k) 


Pi-j-{k)Rj 




ieu 


• If n (k) = fj, + (k), which means that the oriented demand is perfectly balanced among the output 
links, distribute the unassigned demand proportionally to the allocated supply: 

Pi—j (k)Rj 


Pt-j{k +1) = Pt-j{k) + 


1-3 

kj{k + l) 

T±{k+ 1 ) 

fi{k+l) 
• Else, assign: 


E,-« 


feb c _ (fc) 


Pi-j'(k)Rj 




= kj(k) for {*, j,c} 7 ^ {i ,j,c }; 


= 0 ; 


A/3 l c _ i _(fc) = min <( (A:), 


/9» (&) for {*, 4 j- {i ,c }. 

(fc)-R,- E2 =i E-,- ( fc ) 




sf (*) 


£•->(£+ 1) 

^■(fc + i) 
^(fc + i) 


Ay(fc) for {*,j,c} ± {i~,j~,c~}-, 
/3-(fc) for {*, c} yf {*",c"}. 


10. Update sets Uj{k) and V(k): 

Uj(k + 1 ) = Uj(k)\{i- : ^ C _(fc + 1 )= 0 , c=l,...,c}, jGU(fc); 

U(fc + 1) = V(fc)\{j: £/,(£; +1) =0}. 

11. Set k := k + 1 and return to step 2. 


(4.4) 

(4.5) 


(4.6) 

(4.7) 

(4.8) 


2 

m 

1 


4 

$ 

3 


Figure 4: Merge-diverge intersection example. 

Example. Consider the node presented in Figure 4. The node represents a junction between two parallel 
links, or a merge-diverge. Nodes such of these may be used to model locations where drivers can choose to 
shift between two parallel roads or lane groups; the canonical example is an access point for a managed lane 
facility adjacent to a freeway, where drivers may choose to enter or exit the managed lane. In this example, 
we will say that links 1 and 3 represent a freeway’s general purpose (GP) lanes, and links 2 and 4 represent 
a high-occupancy vehicle (HOV) lane, open only to a select group of vehicles. The Low-Occupancy Vehicles 
(LOVs) and HOVs will be modeled as separate commodities, with notations of L for LOVs and H for HOVs. 

Suppose we have input parameters of: 
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Pl4 ~ 0 Pl3 ~ 1 

0E =? P& =? pi = 3 A 

/?! =? P 24 =? P2 = V4 

Sf = 500 5f = 100 R 3 = 600 

S.f = 50 R 4 = 200 

In words, the LOVs are not allowed to enter the HOV lane; HO Vs are allowed, but individual vehicles may or 
may not choose to do so: both links can bring them to their eventual downstream destination so the choice 
of link to take will be made in response to immediate local congestion conditions. 

Let us see how algorithm would assign split ratios for the HO Vs. 

fc = 0 : 


1 . /?f 3 (o)= pm=pm = ~pm= 0 
Pz (0) = P 4 (0) = 1 

A 3 (0) = U 4 {0) = {1,2}; V(0) = {3,4} 

2 . V(0) = {3,4} ^ 0, the algorithm proceeds. 

3. ~sf (0) = 100, sf(0) = 50, sf(0) = 0 

4. Sm = 500; Sf 3 (0) = 5g(0) = S&(0) = S£(0) = 0 

5- 7l3 = 7l4 = V 2 ; 123 = 124 = V 2 ! 7l3 = 1 


Pl3 (0) 


Pi (713 s ! 1 + 7l3 S l) 

sf + sf 


3/4 (1/2 X 100 + 1 X 500) 
100 + 500 


Pu(0) 
P 23 (0) 


3/4 ( 1/2 x 100 + 0 x 500) 


100 + 500 
1/4 (V 2 x 50) 

A) 


= 0.125 


p 24 (0) = 0.125 


0.0625 


0.6875 


6. fi+( 0) 


5 ^( 0 )+sm 

Pl 3 (0)l?3 


Pis( 0) = 0.9848, for all other pairs the ratio equals 0. 

ieu 3 


7. Y (0) = {3,4}; 
j~ =4 


8. 1+4(0) = {1, 2}, as there are no oriented demands to link 4 as of yet. 
(i~ c") = (2 ,H) 

9. n~(0) = 0 


A Pi-j- (0) = A/3^(0) = min 
^(1) =0 + 1 = 1 
Af (1) = 1-0 = 0 


/ 0.9848 x 0.125 x 200 
V 50 x 0.1875 



= min {1, 2.616} = 1 


10. t/ 3 (l) = U 4 {1) = {1}; V(1) = {3,4} 


k= 1 : 
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2. V(l) = {3,4} 

3. sf(l) = 100, sf( 1) = 0, ~Si (1) = 0 

4. S&C0) = 500; Sf 3 (l) = 5g(l) = 5^(1) = 0; S£(l) = 50 

5- 7i3 = 7i4 = V 2 ; 723 = 0; 7 m = 1; 713 = 1 

Pi 3 (l) = 0.6875; p 14 (l) = 0.0625; 

P 23 (l) = 0; p 24 (l) = 0.25 


6 . 


Sh(l)+S&{l) 

Pi 3 (1)-R 3 

g 2 L 4 (l)+^ 4 (l) 

p 24 (l)i?4 


E ft 3 (l) = 0.8333 

ieu 3 

E P<4(1) = 0.3125 

ieUi 


M + (l) 


For the two other pairs the ratio equals 0. 

7. Y(l) = {4}; 
j =4 

8. W 4 (1) = {1} 

{i~c~) = (1, H) 

9. p~( 1) = 0 


A^:._(l)=A^f 4 (l) = min 

/3g(2) = 0 + i/3 = Vs 
-p2 (2) = 1 - 1/3 = 2 /3 


0.8333 x 0.0625 x 200 
100 x 0.3125 



= min{l, i/ 3 } 


V 3 


10. C/ 3 (2) = t/ 4 (2) = {1}; 1/(2) = {3,4} 


k = 2 : 

2. 1/(2) = {3,4} 

3. Sf (2) = 66.67, Sf(2) = 0, sf (2) = 0 

4. S&(2) = 500; ^(2) = 50; S*(2) = 33.33; S»{2) = S£(2) = 0 

5- 7g = Vs; 7u = 2 /s; 7i = 0 ; 7m = 1; 7iV = 1 
P 1 3 ( 2 ) = 2 / 3 ; pu(2) = 1 / 12 ; 

P 23 (2) = 0; p 24 (2) = 0.25 


6 . 


Vf 3 (2)+Vf 3 (2) 

Pi 3 (2)i? 3 

Vf 4 (2)+Vf 4 (2) 

Pl 4 (2)i? 4 


E p ? ; 3 (2) = 0.8333 
ieu 3 

E Pu( 2) = 0.6667 

ieUi 


4 L 3 (2)+Vj(2) 

P 23 ( 2)i? 3 


E P«( 2 ) = 0 

ieu 3 


g 2 V(2)+Vf 4 (2) 

Pm(2)7?4 


E Pm(2) = 0.3333 

ieUi 


M + (2) 
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7. Y( 2 ) = {3}; 

r = 3 

8. W 3 (2) = {1} 

(i~c~) = (1, H) 

9. n~(2) = 0.8333 = /x+(2) 


^ 2 /q v f)00 

«< 3 ) = ° + 2 /3 X 600 + 1 /.2 x 200 X 2 / 3 = ° + ° * X 2 /“ = °- 64 


^ 23 ( 3 ) — V 3 + 


1/12 x 200 


2/3 X 600 + 712 x 200 
/T/ (3) = 2/3 - 0.64 - 0.0267 = 0 

10. £73(3) = £74(3) = 0; 7(3) = 0 


x 2/3 = y 3 + 0.04 x 2/3 = 0.36 


k = 3 : 


2. 7(3) = 0 

The algorithm terminates. The resulting split ratios are (3^ = 0.64, /3|| = 0.36, — 7 /3|| = 0. 


5 Conclusion 


This paper introduced the first order macroscopic traffic model, LNCTM, for modeling multi-commodity 
traffic on road networks. An example of the LNCTM application is the freeway network with special lanes, 
where different vehicle classes have different permissions to enter those special lanes. 

One of the features of the LNCTM is the fundamental diagram in the shape of the “inverse lambda” that 
allows modeling of the capacity drop and a hysteresis behavior of the traffic state in a link that goes from 
free flow to congestion and back (Section 2). 

The cornerstone of the LNCTM is the model of a node with multiple input and multiple output links, which 
computes traffic flows per commodity from input to output links and is governed by the input demand, output 
supply, split ratios that define how incoming flows must be distributed between output links, and input link 
priorities that define the output supply allocation for incoming flows (Sections 3.1-3.2). The proposed node 
model accepts arbitrary input link priorities, and it is shown that the node model presented in Tampere 
et al. ( 2011 ) is the special case of the proposed node model when input link priorities are proportional to 
input link capacities (Sections 3.3). We also analyzed the non-convex nature of the node model, showing its 
difference from the Bliemer node model (Bliemer, 2007) that is formulated as an LP (Section 3.3). 

Next, we discussed the overly restrictive nature of the FIFO rule on the output flows and further generalized 
the node model by relaxing the FIFO condition. Parameters, called mutual restriction coefficients allow 
tuning the node model from no FIFO when they are set to 0 to full FIFO when they are set to 1 (Section 3.4). 
Then, it was shown that this generalized node model produces unique input-output flow allocation that 
maximizes the total output flow of the node, subject to constraints imposed by input demand, output 
supply, split ratios, mutual restriction coefficients and input link priorities (Section 3.5). 

In freeways with activated HOV lane LOVs are confined to the GP lane, whereas HOVs may choose between 
GP and HOV lanes based on traffic conditions, and so split ratios for commodity representing HOV traffic 
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may not be a priori defined, but must be computed at runtime. Therefore, we introduced the algorithm 
for local traffic assignment that computes a priori undefined split ratios at nodes based on input demand, 
output supply and input link priorities (Section 4). 

The above mentioned results are presented in the form of constructive computational algorithms that are 
readily implementable in traffic simulation software. 


Acknowledgements 


We would like to express great appreciation to our colleagues Elena Dorogush and Ajith Muralidharan for 
sharing ideas, Ramtin Pedarsani, Brian Phegley and Pravin Varaiya for their critical reading and their help 
in clarifying some theoretical issues. 

This research was funded by the California Department of Transportation. 


A Notation 


Symbol Definition 
t Timestep index 

T Final timestep 

k Iteration index 

l Generic link index 

C Set of all links 

v Generic node index 

A f Set of all nodes 

i Index of links entering a node 

M Number of links entering a node 

j Index of links exiting a node 

N Number of links exiting a node 

c Vehicle commodity index 

C Number of vehicle commodities 

rii Density of vehicle commodity c in link l 

ffj Flow of vehicle commodity c from link i to link j 

Fi Capacity of link l 

vi Speed of flow in link l 

vf Freeflow speed of link l 

u>i Congestion wave speed of link l 

nj Jam density of link l 

nj Low critical density of link l 

High critical density of link l 
f3fj Split ratio of c vehicles from link i to link j 

i)jj. Mutual restriction coefficient of link j onto link j 1 for link i 

61 Congestion metastate of link l 

Si Demand for commodity c of link l 

Si(k ) Adjusted demand for commodity c of link l as of iteration k 

Ri Supply of link l 

Pi Priority of input link i 

U(k ) Set of input links whose flows have yet to be fully determined as of iteration 
k 


Used in section(s) 
2 
2 

3,4 

2 

2 

2 

2 

3,4 

3,4 

3,4 

3.4 

2.3.4 
2,3,4 
2 

2,3 

2,3 

2 

2 

2 

2 

2 

2 

2,3,4 

2,3 

2 

2,3 

3 

2.3 

3.4 
3 
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V(k ) Set of output links whose flows have yet to be fully determined as of iteration 3 

k 

pi(k ) Adjusted priority of link i at iteration k 3,4 

Rj(k) Adjusted supply of link j at iteration k 3 

Uj(k ) Set of input links contributing to link j whose flows are undetermined as 3 

of iteration k 

Pij(k) Oriented priority from link i to j at iteration k 3,4 

a,j(k) Restriction term of link j at iteration k 3 

dj * (fc) Smallest (most restrictive) restriction term at iteration k 3 

ctj(k) Reduction factor of link j at iteration k 3 

U(k) Set of input links whose demand can be fully met by downstream links at 3 

iteration k 

B Set of commodity movement triples (i, j , c) whose split ratios are known 4 

B Set of commodity movement triples (i, j, c ) whose split ratios are unknown 4 

/3° Unassigned portion of the split ratios of commodity c from link i 4 

Unassigned demand of commodity c from link i 4 

V(k) Set of output links to which unassigned demand may still be assigned as of 4 

iteration k 

p + (k) Largest oriented demand-supply ratio at iteration k 4 

pT (k) Smallest oriented demand-supply ratio at iteration k 4 
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